% Copyright (C) 2009,2010,2011,2012  Marco Restelli
%
% This file is part of:
%   LDGH -- Local Hybridizable Discontinuous Galerkin toolkit
%
% LDGH is free software: you can redistribute it and/or modify it
% under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% LDGH is distributed in the hope that it will be useful, but WITHOUT
% ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
% or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public
% License for more details.
%
% You should have received a copy of the GNU General Public License
% along with LDGH. If not, see <http://www.gnu.org/licenses/>.
%
% author: Marco Restelli                   <marco.restelli@gmail.com>


function [p,e,t] = new1d_grid(ne,a,b)
% [p,e,t] = new1d_grid(ne,a,b)
%
% Define a 1D grid with ne elements on the interval [a,b].

 if(b<a)
   x = b;
   b = a;
   a = x;
 end

 p = linspace(a,b,ne+1);
 t(1,:) = [1:ne];
 t(2,:) = [2:ne+1];
 e(:,1) = [ne+1;0;0;1];
 e(:,2) = [   1;0;0;2];

return

